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Abstract We derive the order parameter of the chiral Potts model, using the method of Jimbo 
et al. The result agrees with previous conjectures. 

PACS numbers: 05.50.+q, 64.60.Cn, 75.10.Hk 

The solvable chiral Potts model is one of iV-state spins with nearest-neighbour 
interactions on a planar lattice C P It is like other solvable models in that its inter- 
actions satisfy the star-triangle relations. However, it is unlike most of them in that 
it does not have the "rapidity difference" property. This makes the model mathe- 
matically much more difficult. The free energy of the model has been obtained using 
functional transfer matrix relations, [21 121 E] but the order parameters (spontaneous 
magnetizations) have so far defied calculation. In particular, the corner transfer 
matrix method appears to fail completely. |Hj 

This is despite the fact that there is an elegant and eminently believable conjec- 
ture for these order parameters. 13 E] 

Here we derive the order parameters and verify the conjecture. The method is 
quite similar to the "inversion relation" method for calculating the free energy. [HI 
Both involve assuming certain analyticity properties of a ^^T2{tq)" model (which is 
closely related to the superintegrable chiral Potts model). 

We take C to be the square lattice, drawn diagonally as in Fig. [T] Each spin i 
takes one of possible states, labelled 0, . . . , A^ — 1. It interacts with a neighbouring 
site j with Boltzmann weight Wyh{cri — aj) or W^hicTi — aj) if i is below j and the 
edge is in the SW-NE or SE-NW direction, respectively. 

Here v, h each denote a set of rapidity variables. Let p = {xp,yp, fip,tp}, where 
Xp, Up, fip, tp are related to one another by 

+ < = Ml + x^y^!) , kx^ = l- k'/ix^ , ky^ = 1 - k'ft^ , (1) 



and 



= {i-k'y/^ , o<k,k'<i . (2) 
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Here k, k' are constants, the same at all sites of the lattice, and k' is a "temperature- 
like" variable, being small at low temperatures. Here we take < /c' < 1, which 
means that the system is ferromagnetically ordered. It becomes critical when /c' — > 1. 

Rapidities such as p, g, f , h are associated with the dotted lines of Fig. [T] and 
may vary from line to line. If the edges («, j), (i, k) of C are intersected by rapidity 
lines f , /i, ordered as in Fig. HI then the edge weight functions are Wvhicn — o'j), 
Wyh{ai - ak) where 

Wp,{n) = n , (3) 

j = l Up ~ ^•'■I^q 

Wp,{n) = n ^^^^^^^ , (4) 

where u = exp{27ii/N) and v = {xy,yy, fiy,ty}, h = {xh,yh, ^h,th} are rapidity 
sets satisfying (P). 

If Xp, Xq, Hp, Hq, uJXp all llc ou the unit circle in an anti-clockwise sequence, then 
Wpqiji), Wpqin) are real and positive. We refer to this as the physical case. 
Hold the spin a in Fig. Q fixed. Then the partition function is 

Z{a) = Y.X{WyS-j)X{WyHU~k) , (5) 

cr 

the products being over all edges of each type (with their appropriate rapidity vari- 
ables), and the sum over all values of the other spins, the boundary spins being set 
to zero. The unrestricted partition function is Z = Z{0) -|- ■ ■ ■ -|- Z{N — 1). 

We have followed Jimbo et al ^0] and cut the horizontal rapidity line immediately 
below a, giving the left (right) half-line a rapidity p (q). Define 

Fpq{a) = Z{a)lZ . (6) 

We expect this ratio to tend to a limit when the lattice is large. It will not depend 
on the "background" rapidities v^h because of the star-triangle relation, [T] which 
allows us to move any of the w , h rapidity lines infinitely far away from the spin a. 

However, we cannot move the half- lines p, g away from a, so Fpq[a) will indeed 
be a function of p, q. It is the probability that the central spin has value a. An 
exceptional case is when q = p., when the cut disappears and the recombined line 
can be moved to infinity, so 

Fpp{a) = independent of p . (7) 

What we can do is rotate p, q round a and then interchange them, which gives 
us some functional relations satisfied by Fpg(a).||I^ If we define 

N-l 

Fpq{r) = J2 ^'^'^ Fpqia) , (8) 

a=0 

Gpq{r) = Fpq{r)/Fpq{r - 1) , (9) 
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Figure 1: The square lattice with the cut rapidity hue below a. 



then for the infinite lattice we obtain the relations 



N 

I[Gpq{r) = 1 



G 



R-^q,Rp 



(r) , 



r=l 



GMp,q{r) = Gp^M-^qir) = Gpg{r + 1) . 
Here R, M are automorphisms that act on the rapidities: 

{XRp,yRp,fiRp,tlip} = {yp,UJXp,l/ flp,Ujtp} 



(10a) 
(10b) 

(10c) 

(lOd) 
(lOe) 

(11) 



{xMp,yMp, fJ'Mp,tMp} — {Xp,yp,UJfip,tp} . (12) 

We regard tp as an independent complex variable and Xp,yp,fip as determined 
from it by (0). They are multi- valued functions of tp-. to make them single- valued 
we must cut the tp-plane as in Fig. |21 There are cuts Cq, . . . , C^^i, where Cj lies 
on the radial line arg(tp) = 27Tj/N. 
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Figure 2: The tp-plane for = 3, showing the cuts Cq, Ci, C2. 

The case we shall be interested in is when |/ip| > 1 and arg(xp) is between 
—7r/{2N) and 7c/{2N). Then Xp lies in a small approximately circular region TZq 
round Xp = 1, while i/p lies in a plane with corresponding approximately circular 
holes TZq, . . . , T^AT-i surrounding the points 1, a;, ... , We say that p lies in the 

domain T>i. 

In the low temperatures limit k' ^ and Cj, TZj both shrink to the point 
Then Xp — 1 and yp can lie anywhere except at a root of unity. We take p,q,h e T>i 
and = 0(1), Xjj — Vv — 1- Then (n) , VT^/i (n) are small unless n = (mod 
A^), which is the usual low-temperature case. The sum in ^ is dominated by the 
contribution from all spins other than a being zero, and we obtain 

_ k'\fip/f,,r f.i-u^ 

Let us take q to be related to p by 

Xq = Xp , yg = uyp , fig = fip . (14) 

Then ()13p simplifies to 
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Thus to leading order in k' we see that Fpq{a) is an analytic and bounded function 
of tp except at tp = 1. This is consistent with Fpg{a) and Gpq{r) having branch cuts 
(for non-zero k') on Cq, but there is no indication of branch cuts at Ci, . . . , Cn-i- 

We argue that this is exactly true for sufficiently small non-zero values of k' . If 
we rotate and reverse the left half-line p anticlockwise below a, as in ^T], it becomes 
a left-pointing line with rapidity p' = R~^p} (To ensure that at low temperatures 
the dominant contribution continues to come from all non-central spins being zero, 
we should reduce /i^ to order fc', to order unity, and take Xh — Uh — 1-) 

The half-line p' now lies immediately below g, and from (|THl 

Xq = yp' , yq = Uj'^Xp' , /ig = . (16) 

This is precisely the condition for the combined weights of the two half-rows p', q 
(summed over intervening spins) to be those of the T2{tpi) model. (Eqns. 2.38 to 
3.48 of 0, with k = 0, i = 2.). For a finite lattice, if there are s spins to the left of 
a, then Z{a) depends on p only via tp, and is a polynomial of degree s in tp. 

Further, when k' is small and a = 0, this polynomial is (1 — tp)^, and is small (of 
order fc'^) when < a < A^. Hence Gpq{r) is the ratio of two polynomials, each of 
degree s and tending to (1 — tpY as k' — > 0. 

By continuity, for sufficiently small but non-zero k' the zeros of the polynomials 
must be close to tp = 1, i.e. tp/ = l/uj. 

This is the same behaviour as the free energy of the T2{tpi) model. This leads us 
assume, corresponding to Assumption 2 of P, that Gpq{r) is analytic in the cut tp 
plane of Fig. |2l except only for the branch cut Cq? 

Now look at (|lUc|) . Using (jl4j) it becomes, for r = 1, . . . , — 1, 

Xp^Gpqij) = x=^Gp^q{r) , (17) 

where p = R^^q,q = Rp. It follows that p,q satisfy the relation (fT^ and Xp = 
yp, yp = Xp, fip = l//ip. The LHS of this equation is therefore the same as the RHS, 
with the same value of tp, but with Xp, yp interchanged and fip inverted. 

This is what happens if one crosses the branch cut Co in Fig. |21and then returns 
to the original value of tp. Thus (fT7j) is equivalent to the statement that Xp^Gpq{r) 
does not have a branch cut at Cq. 

Write Gpq{r) as Gr(tp) and consider the function 

Lritp) = X;'Gr{tp)Gr{ujtp)---Gr{uj''-Hp) . (18) 

For r 7^ the factor Xp^Gr{tp) has no cut on Cq. From our assumption, neither 
do any of the other G factors. Hence Lritp) has no such cut. It is unchanged 
by tp uotp, so has no cuts at any of the Cj. When tp,yp,fip oo (their ratios 
remaining finite and non-zero) the Boltzmann weights remain finite, so we expect 

^We can then interchange p' with q, which gives the symmetry (I10b|l . 

^One can obtain additional evidence for the analyticity near Cq from the figure obtained by 
rotating p clockwise to become a half-line Rp lying above ».[T^ 
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Lr{tp) to be bounded at infinity. From Liouville's theorem it is therefore a constant, 
so 

Lr{tp) = Cr , 0<r <N , (19) 

Cr being a constant.^ 

When tp, Up = 0, then q = p and Xp = k. If tp, yp, fip — > oo, then Xp = 1/k and 
q = M~^p. From dZj) and (fiTlejl it follows that 

Cr = k-'/^'Gppirf = fci/^ Gpp{r + if , (20) 

for r = l,...,A^ — 1. 

Eliminating Cr and using ()10d|l . we obtain 

Gpp{r) = (21) 

for r = 1, . . . , N. 

The order parameter of the chiral Potts model is 

(^™) = Fp,ir)/Fp,iQ) 

= Gpp{l)---Cpp{r) , (22) 

so we have 

(CU™) = k'-^N-r)/N^ (23) 

for r = 0, . . . ,N. This is the result previously conjectured on the basis of series 
expansions (eqn. 1.20 of Ref. jH]). We expect all functions to be analytic in the 
physical ferromagnetically ordered regime (with positive real Boltzmann weights) so 
our working and results should remain true throughout < fc' < 1. The magnetic 
critical exponents are jSr = r{N — r)/{2N'^). 

In we go on to obtain Gpq{r) from a Wiener-Hopf factorization of Xp. The 
result is a special case of the T2{p') free energy. For r = A^, it and Lr{tp) can be 
obtained from ()10d|) . For = 3 we have verified that these results are consistent 
with previously obtained series expansions (eqns. 48 - 52 of Ref. |E1)- 
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